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Abstract. Lya forest lines in QSO spectra have typical widths of 20 — 50 km s~^. 
Low column density absorbers in cosmological simulations are large, diffuse struc- 
tures, and the Hubble flow across the spatially extended absorber is usually the dom- 
inant contribution to the width of its associated absorption line. Thermal broadening 
is unimportant over most of the spectrum, and peculiar velocities tend to make ab- 
sorption features narrower rather than broader. As a consequence of Hubble flow 
broadening, there is a close relation between local Lya optical depth and local neu- 
tral hydrogen density, which is well approximated by the Gunn-Peterson formula. 
The physics that governs the unshocked intergalactic medium leads to a tight corre- 
lation between the neutral hydrogen density and the underlying gas and dark matter 
overdensity. For many purposes, it is simpler to regard a Lya forest spectrum as a 
continuous, non-linear map of the density field rather than a collection of discrete 
lines. This continuous field view of the Lya forest can be applied to measurement of 
the baryon density parameter, testing of cosmological models, and robust determi- 
nation of the shape and amplitude of the primordial mass power spectrum. 

In the last few years, the study of the Lya forest has undergone several 
observational revolutions: the extension to low redshift via HST, the probe 
of internal structure from spectra along neighboring lines of sight, the ex- 
traordinary detail provided by Keck HIRES data, and the clear detections 
of metals associated with low column density HI absorbers. The field has 
also undergone a theoretical revolution, driven by hydrodynamic cosmological 
simulations. These allow one to predict properties of the Lya forest from a pri- 
ori theoretical models motivated by independent considerations of large scale 
structure, the cosmic microwave background, and galaxy formation. Since the 
pioneering numerical study of Cen et al. there have been more than 30 
papers using cosmological simulations to investigate QSO absorption phenom- 
ena, by several independent groups. The resulting picture of the Lya forest 
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has features in common with some earUer models, especially the fluctuating 
intergalactic medium (IGM) scenario of Bi |Q . 

One distinctive feature of this cosmological picture of the Lya forest is 
the low density of the absorbing structures. Typical marginally saturated 
lines (A^Hi ^ 10^"' cm^^) arise in gas whose density is a few times the cosmic 
mean or less. Weak lines (iVni ^ 10^'^ cm~^) often occur at local maxima 
that lie below the global mean density. This low density has a number of 
important consequences. Most absorption arises in structures that are still 
expanding with residual Hubble flow. These absorbing systems are usually 
far from dynamical, hydrostatic, or thermal equilibrium. The low density 
implies a low recombination rate and thus a low neutral fraction (typically 
~ 10~^ — 10"*), so the neutral hydrogen revealed by the observed Lya opacity 
is only the tip of a much larger iceberg. 

The low density of the absorbing gas also means that absorbers must be 
physically large in order to produce the observed column densities. The large 
size implies that the Hubble flow across an absorber can be substantial. Indeed, 
for a typical low column density line in the simulations, the primary contribu- 
tion to the line width (6-parameter) comes from Hubble flow. This situation 
contrasts with that in traditional conceptions of the forest, where lines are 
assumed to be broadened by thermal motions of the gas or by "turbulent" 
motions of cloudlets. Gravitationally induced peculiar velocities do distort the 
lines in cosmological simulations, but these coherent flows are not at all like 
Gaussian turbulence, where there is a large dispersion in the velocities at a 
given spatial position. Furthermore, because most lines with A^hi ^ 10^^ cm~^ 
arise in moderately overdense regions that are expanding slower than the Hub- 
ble rate, peculiar velocities on average have the effect of narrowing Lya forest 
lines, not broadening them. 

Figure 1 illustrates these points using spectra extracted along twelve ran- 
dom lines of sight through a hydrodynamic simulation of the "standard" CDM 
model (SCDM, with n = 1, h = 0.5, erg = 0.7) at z = 3 (see § and [| for 
details). Solid lines show the full Lya absorption spectra. Dotted lines show 
the spectra with no thermal broadening — they are computed by artificially 
setting the gas temperature to zero (without changing the neutral fraction). 
There are a handful of sharp features in the non-thermally-broadened spectra 
that are smoothed away in the full spectra. However, in most regions the 
dotted and solid lines are barely distinguishable, demonstrating that thermal 
broadening usually does not contribute significantly to the width of the ab- 
sorption features. The dashed lines in Figure 1 show spectra with thermal 
broadening but no peculiar motions — they are extracted along the same lines 
of sight after setting the peculiar velocities of all gas particles to zero. Compar- 
ing the solid lines to the dashed lines shows that peculiar velocities do shift the 
positions and distort the shapes of individual absorption features, but they do 
not make them systematically broader. Indeed, as expected from the physical 
argument above, the lines in the full spectra tend to be somewhat narrower 
than the lines in the spectra with no peculiar velocities. 




Figure 1: Solid lines show spectra along twelve randomly chosen lines of sight 
through a simulation of the standard CDM model, at z = 3. Dotted and 
dashed lines show spectra along the same lines of sight with no thermal broad- 
ening and no peculiar velocities, respectively. 

Studies of QSO pairs [0, ||, 0] provide observational evidence for Hubble 
flow broadening of Lya forest lines independently of cosmological simulations. 
The inferred transverse coherence scale of the absorbers. It ~ 150/i~^ kpc at 
z ~ 2, corresponds (for = 1) to a line-of-sight extent H{z)lt ~ 80 km s~^, 
considerably larger than the ^ 25 km s"^ 6-parameter of a typical forest line 
[ p^ . Nonspherical absorbers would be preferentially intercepted "face on," but 
unless the absorbing structures are highly flattened the observed transverse 
scale implies that Hubble flow across them must make a major contribution to 
the line width. The transverse coherence could in principle be a signature of 
clustering of small clouds rather than a physical scale of large absorbers, but 
the nearly perfect coincidence of lines towards small-separation gravitational 
lens pairs |l^ argues against the clustering interpretation, and the detailed 
match of absorption features shown by Ranch in these proceedings seems to 
rule it out definitively. 

Many papers have remarked on the low density of the absorbing gas in 
cosmological simulations (e.g., [|[ ^ ||, ^ ^). The issue of Hubble flow 
broadening has received less attention, but its implications are perhaps even 
more profound. For a start, it means that the profile of a Lya forest line 



shows a line-of-sight density profile through the absorbing structure, albeit 
one that is non-linear and distorted by peculiar motions. This is not the case 
in the thermal broadening picture, where the absorber itself is compact and the 
wings of the line arise from high velocity atoms. In the cosmological picture, 
line wings show the absorbing structure itself fading into the background, like 
mountains into foothills. 

Another consequence of Hubble flow broadening is that the Gunn-Pctcrson 
[D formula. 



provides a good approximation to the relation between local Lya optical depth 
and the local space density of neutral hydrogen. In a thermally broadened, 
compact cloud model, by contrast, the optical depth is lower than tqp at the 
line center (the redshift space location of the dense cloud) and higher than tqp 
in the line wings (where there is no gas at the corresponding redshift space 
position). Cosmological simulations imply that the Lya forest can be viewed 
as a fluctuating Gunn-Peterson effect, produced by an inhomogeneous, diffuse 
intergalactic medium. 

What turns the fluctuating Gunn-Peterson idea from a novelty into a pow- 
erful conceptual tool is the simplicity of the physics that governs the ionization 
state of the low density gas. This gas is in photoionization equilibrium, so 
the neutral hydrogen density is nni oc p^T~^ '^ /T, where F is the photoion- 
ization rate and the T~^''^ factor accounts for the temperature dependence 
of the hydrogen recombination coefficient near T - 10-* K. The interplay 
between photoionization heating and adiabatic cooling leads to a tight rela- 
tion between temperature and density, which can be well approximated by 
a power law, T = Toip/py [|, 0. The values of Tq and 7 depend on the 
UV background spectrum and reionization history and can be computed semi- 
analytically typically Tq - 6000 K and 7 ~ 0.3 - 0.6. With this physical 
reasoning, equation (^ can be converted to a formula we can describe as the 
fluctuating Gunn-Peterson approximation, 



where (3 = 2 — O.77, p/p is the overdensity at the position where the redshift 
(cosmological plus peculiar velocity) is Aobs/Aa — 1, and dVios/dx is the deriva- 
tive of the line-of-sight peculiar velocity at the same position. The peculiar 
velocity term accounts for the mapping from real space to redshift space. In 
principle, p/p here refers to the gas overdensity, but because the temperature 
is low, pressure gradients are small compared to gravitational forces, and the 
gas traces the dark matter quite well. 

Equation (|^) is valid if all gas lies on the temperature-density relation 
and thermal broadening and coUisional ionization can be ignored. The ap- 
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proximation breaks down when pj'p ^ 10, but these regions occupy a small 
fraction of the spectrum. They are responsible for high column density lines, 
and the general description presented above begins to break down for TVhi ^ 
10^^ — 10^^ cm~^. This description might also become less accurate at low 
redshifts; we have not examined the Lya forest at z < 2 with our simulations. 
We should also note that some low column density lines at high z arise in 
shock heated gas and are thermally broadened. 

A second useful approximation arises from ignoring peculiar velocities, set- 
ting dV\os,ldx = in equation (H), so that there is a one-to-one relation be- 
tween optical depth and overdensity. Figure 6 of j|] shows that the t — p 
relation remains tight in simulated spectra even when peculiar velocities, ther- 
mal broadening, shock heating, and coUisional ionization are all taken into 
account. The distribution of Lya optical depths P(t) is directly observable 
from high resolution QSO spectra and one can use this second approxi- 
mation to write the mean density of the "warm" IGM that produces the Lya 
forest in terms of an integral over this distribution, PwiGM = /q p{T)P{T)dT. 
After some manipulation, one obtains the density parameter of this warm IGM 
component jl^, 17 , 
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For the fiducial value of the optical depth integral, we have used a value inferred 
from the observations of |^3| at z = 3. The implied f^wiGM is a substantial 
fraction of the baryon density parameter f2f, allowed by big bang nucleosyn- 
thesis, indicating that most of the baryons in the universe at z ~ 3 resided in 
the Lya forest, as the simulations predict. 

The relation between the underlying mass distribution and the number 
density of Lya forest lines in a spectrum may be physically complex, and 
it is sensitive to the details of the observational procedures and the method 
used deblend absorption features. However, the fluctuating Gunn-Peterson 
approximation implies that the relation between mass density and observed 
flux is direct and simple. If one wants to use the Lya forest to test theories 
of structure formation, it is best to abandon lines altogether and treat the full 
observed spectrum as a continuous field. Statistical properties of this "flux 
field" are directly related to statistical properties of the underlying density 
and velocity fields, which are basic predictions of cosmological models. We 
are currently studying a variety of statistical measures similar to those used in 
large scale structure analyses, applying them to simulations in order to assess 
their sensitivity to different properties of the initial fiuctuations and to values 
of cosmological parameters. We have also developed and tested a method to 
recover the shape and amplitude of the primordial mass power spectrum P{k) 
from Lya forest data |5|, again motivated by the "continuous field" point of 



view. The 3-d flux power spectrum has the same shape as the mass power 
spectrum on large scales, and the normalization can be determined by evolv- 
ing numerical simulations with this initial P{k) shape until they reproduce 
the observed power spectrum of the QSO flux. Imposing the observed mean 
Lya opacity as a constraint makes the derived P{k) normalization insensi- 
tive to the choice of cosmological parameters, ionizing background spectrum, 
or reionization history. This approach thus neatly circumvents the uncertain 
physics of galaxy formation and "biasing," which complicates the interpreta- 
tion of power spectra measured from galaxy redshift surveys. Application to 
existing samples of QSO spectra should soon yield the power spectrum of mass 
fluctuations in the high redshift universe. 
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